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In statistical mechanics, it is well known that finite-state classical lattice models can be recast as 
quantum models, with distinct classical configurations identified with orthogonal basis states. This 
mapping makes classical statistical mechanics on a lattice a special case of quantum statistical me- 
chanics, and classical combinatorial entropy a special case of quantum entropy. 

In a similar manner, finite-state classical dynamics can be recast as finite- energy quantum dynam- 
ics. This mapping translates continuous quantities, concepts and machinery of quantum mechanics 
into a simplified finite-state context in which they have a purely classical and combinatorial inter- 
pretation. For example, in this mapping quantum average energy becomes the classical update rate. 

Interpolation theory and communication theory help explain the truce achieved here between perfect 
classical determinism and quantum uncertainty, and between discrete and continuous dynamics. 



INTRODUCTION 

In this paper we discuss a mapping between classical 
and quantum systems that lets us regard quantum dy- 
namics as a generalization of finite state classical dynam- 
ics, and that allows us to identify equivalent quantities 
and concepts in classical and quantum systems. 

A similar mapping has long been known in statistical 
mechanics [T] that establishes classical lattice models and 
their combinatorial entropy as simple examples of quan- 
tum statistical mechanics. 

There is an obvious candidate for the comparable dy- 
namical mapping: classical computations are equivalent 
to a subset of quantum computations [5J- Most work 
on quantum computation is, however, based on hybrid 
classical/quantum models in which macroscopic classical 
operations control the sequencing of quantum operations. 
Such systems do not provide a purely quantum target for 
a classical/quantum mapping. Instead, early work show- 
ing that autonomous quantum systems can perform clas- 
sical computation [3] forms the basis for the dynamical 
mapping presented here. 

This mapping allows physical quantities such as energy 
and momentum to be identified with finite-state classical 
quantities, with the aid of classical interpolation theory. 
Related issues are addressed in [3] , but a general dynam- 
ical mapping is not provided there. 

As a preliminary to discussing dynamics we first review 
a canonical method for mapping classical lattice models 
onto quantum lattice models in statistical mechanics. 

STATISTICAL MECHANICS 

In statistical mechanics, it is well known that classical 
lattice models can be recast as quantum models, with 
distinct classical configurations identified with orthogo- 
nal basis states PQ. 

Consider, for example, the well known ferromagnetic 
2D Ising model. In this model each of M lattice sites in a 



square lattice is occupied by a classical two-state "spin," 
and each state S n of the N = 2 M possible configurations 
of the lattice is assigned a classical configurational en- 
ergy ^classical fagfc depends only on how many pairs of 
adjacent lattice sites have the same spin value and how 
many have opposite values. 

A quantum lattice model corresponding to such a clas- 
sical lattice model can be constructed by identifying each 
of the N distinct classical states S n with a distinct basis 
vector \n) in an N dimensional Hilbert space. A hamilto- 
nian operator H is defined by taking each configuration 
state \n) to be an energy eigenstate of H with energy 
eigenvalue ^classical . 

H \n) = ££ lassical |„) . (1) 

In quantum statistical mechanics the energy eigen- 
states are also eigenstates of the density operator p, with 
eigenvalues that give the statistical weight to attach to 
each energy eigenstate. For example, for a canonical en- 
semble of quantum mechanical systems, p is proportional 
to e~^ H . From this becomes the usual classical Boltz- 
mann factor when applied to a configuration state \n), 
and quantum statistical mechanics reduces to classical. 

CLASSICAL DYNAMICS 

Since the definition of H used in the statistical mechan- 
ics mapping makes each classical configuration a time- 
invariant state under unitary time evolution, we use a 
different definition of H to emulate classical dynamics. 

Finite-state dynamics 

An invertible classical finite-state dynamics is a dis- 
crete idealization of classical dynamics [5] . Perfect digi- 
tal degrees of freedom are updated at discrete times ac- 
cording to a sequence of invertible transformations. The 
total amount of state in the system, including that used 
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to define the dynamics, is finite. Here we take the time 
between update events to always be r, so that the system 
is updated at the constant rate v = 1/r. 

The finite set of possible configurations of the system 
is partitioned by the invertible dynamics into a collec- 
tion of disjoint dynamical orbits, with each dynamical 
orbit consisting of a set of configurations that turn into 
each other under the dynamics (c/. [SIED- For each dy- 
namical orbit d the number of configurations in the 
orbit determines the period = tN^ of the orbit. One 
configuration of each orbit is labeled with the integer 0. 
The configuration obtained from by one update step is 
labeled 1, and so on. 

We identify configuration n of dynamical orbit d with 
the basis state \n,d). Because the orbit is periodic, 
\0,d) = \Nd,d). This mapping identifies each possible 
configuration of the classical dynamics with a basis state: 
we call this the configuration basis. 



Hamiltonian dynamics 

Given an invertible classical finite-state dynamics, we 
construct a continuous quantum hamiltonian dynamics 
isomorphic to the classical dynamics at regularly-spaced 
times. We begin by defining a discrete Fourier trans- 
formed set of basis states. Let 
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for integer m, where 11 E " is the name of the new basis. 
The inverse transformation is 
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\n,d) = -±= J2 e- 2 ™/^|i?:m,d) 
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We define a hamiltonian H by assigning the \E: m,d) 
states to be its energy eigenstates and E m ^ = m h/Td to 
be the corresponding energy eigenvalues |26j : 



H \E : m, d) = m — \E : m,d) . 
Td 



(4) 



If we let U = e lHT / h be the time evolution operator for 
the time interval r, then 
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Average energy 



The configuration state ^ is a uniform superposition 
of all Nd energy eigenstates \E: m,d) with eigenvalues 



mh/Td, and so the average energy is 



E = 



h{N d - 1) 
2T d 



(0) 



We've taken E d = in the construction above, but the 
fact that the system has a harmonic-oscillator-like en- 
ergy spectrum suggests that we should really add h/2Td 
to all the energy eigenvalues. This is in fact the small- 
est energy allowed by quantum distinguishability bounds, 
assuming the ground state energy of a much larger sys- 
tem encompassing this one sets the zero of the energy 
scale Adding h/2Td makes the average energy ([6| 
independent of Td, 



E = 



hv 
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This is the least possible average energy compatible with 
a dynamics that traverses distinct states at the average 
rate v [9]. Thus our construction is energetically ideal, 
and the average energy is identified with the classical 
update rate of the finite-state dynamics. 

If a lattice dynamics is updated sequentially — one lo- 
cation at a time in a repeating cycle — the frequency with 
which a given location is updated determines a local en- 
ergy. Total update frequency (total energy) is the sum 
of the local frequencies [27]. Different kinds of updates 
{e.g., ones involving particle or bond motion, and ones 
that don't [H1H3]) define different kinds of energy 0]. 

For a large system with a very long period, h/2Td ~ 0, 
and so for simplicity we will revert to taking Eq ^ = in 
the remainder of the discussion. 



BANDLIMITED STATES 

We have provided a prescription for constructing a 
continuous-time quantum hamiltonian description for 
any invertible classical finite-state dynamics — turning 
discrete-time models into continuous-time models. This 
construction can be regarded as an application of band- 
limited interpolation theory [14] . 



Bandlimited dynamics 

Let us choose our unit of time such that r = 1, so that 
our configuration basis states are simply the states seen 
in the dynamics at integer values of time starting from 
\0,d). At a continuous moment of time t the state is 



\t,d) =e~ lHt/h \0,d) 
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multiples of 1/T, and if the lowest frequency is 0, then 



FIG. 1: \S{N, u)\ probability distribution (solid) versus nor- 
malized gaussian of unit height (dashed). 



which is just ^ with t replacing n. We can express 
the continuous-time state \t, d) as a function of the Nd 
integer-time states \n,d) by replacing \E: m,d) with its 
definition 
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where 
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S(N, ll)=Ij] e 2^mu/N 



(10) 



The function S(Nd, n — t) equals the Kronecker delta 5„ it 
for integer values of t between and Nd — 1 but is also 
defined for non-integer values. S(N, u) is a periodic ver- 
sion of the sine function |15j , which is the foundation of 
bandlimitcd interpolation theory: S(N, u) = 1 for integer 
values of u that equal modulo N and S(N, u) = for 
other integer values of u. In fact, if we sum the geometric 
series we recover sine times a phase for large N, 



N 



lim S(N, u) = e* 



(11) 



A portion of the probability distribution ^(JV, u)| is 
shown in Figure [I] for TV = 100 (solid). Near its cen- 
ter it is approximately gaussian (dashed). 



Reconstruction from samples 

Using S(N,u), any periodic function f(t) with period 
T and a bandlimited Fourier spectrum with N frequen- 
cies can be reconstructed from N equally spaced samples. 
Because of the periodicity all frequencies must be integer 



m= 



27rimt/T 
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for some set of a m . Using t = T/N, f(t) is also given by 
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fit) = /( nr ) 
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(13) 



This is obviously true at the N sample times t — nr 
and so it must be true at all times, since S(N, — — n) 
is composed of the same frequency components as f(t), 
and the N coefficients a m are completely determined by 
the values of f(t) at the N sample times t = nr (in fact, 
the the Fourier transform of the /(nr)'s). If the 

lowest frequency is k/T rather than zero, use Sk(N, u) — 
e 2^ku/N S ( N ^ u ) i ns t ea d of S(N, u) above. 

Thus ([9]) can be regarded as an exact reconstruction of 
a continuous but bandlimited dynamics in Hilbert space 
from Nd samples. The bandlimit on the energy spec- 
trum erases the distinction between continuous-time and 
discrete-time dynamics (and field operators [16]). since a 
bandlimited periodic function is completely determined 
by a finite number of sample points. 

If g{t) has the same period and bandwidth as fit) (per- 
haps with a different lowest frequency) then ( 13 ) implies 
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dtfit)git) = ±J2finr)ginT), 



(14) 



and so a bandlimit also erases some of the distinction be- 
tween continuous and discrete analysis of the dynamics. 

CONTINUOUS ISOMORPHISM 

Rather than just have integer time states of a classical 
finite-state dynamics correspond to integer time states of 
a quantum finite-energy dynamics, we can also extend 
the classical finite-state dynamics to intermediate times 
and have the two systems be isomorphic at all times. 

Continuously extended dynamics 

In classical finite-state lattice dynamics it is often use- 
ful to imagine that, when a 1 representing a particle hops 
from one lattice site to another, it moves continuously 
in between. This extension of the dynamics allows us 
to extend classical-mechanical conservations associated 
with continuous spatial symmetries to discrete particle 
motion in order to define, for example, momentum con- 
serving lattice gases [2"5] , 

Continuously extended lattice dynamics have a con- 
tinuous evolution in both time and space but, at every 
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moment, only a finite amount of state: if there are n 
spots in space that can have a 1 or not at integer times, 
there are still only n spots that can at non-integer times. 
Since these n bits don't change their values while they're 
moving between integer locations, the non-integer-time 
states are really just a fixed sequence of rearrangements 
of the bits of the integer-time state. These extra inter- 
mediate states are distinct classically since the bits are in 
different places but they are redundant informationally. 

Note that a continuously extended lattice dynamics 
can still be described as a repeated cycle of local updates, 
but in this case each update moves a bit only infinitesi- 
mally. After any finite interval of time all of the bits will 
have moved by equivalent amounts. 



Continuously extended isomorphism 

In a continuously extended classical lattice dynamics, 
any unit-time separated sequence of states provides a 
complete description of the logical dynamics: since the 
bits of state don't change between integer times, exactly 
when we sample them doesn't matter. 

Similarly, any unit-time separated sequence of states 
from the continuous unitary evolution ([8| constitute a 
complete orthonormal basis set, since ^ implies 



(t',d\t,d) =S(N d ,t'-t) 



(15) 



Thus we are free to define a distinguished basis at any 
time t consisting of the unit-time separated set of Nd 
states from the evolution ^ that includes the current 
state \t,d). If we identify these basis states with cor- 
responding unit-time separated classical configurations, 
then the classical and quantum dynamics are isomorphic 
at all times. 

In analyzing finite-state dynamics, the \t, d)'s act much 
like a complete continuous basis since, again from Q, 

dt \t,d)(t,d\ = \n,d)(n,d\ = I . (16) 

71=0 



Moreover, the inner product (151 acts like a Dirac delta 



function in an integral with a bandlimited function f(t). 
From full, 



dtf(t) (t',d\t,d)=f(t r ) 



(17) 



The continuously extended isomorphism can be used 
to compute average values for operators, such as momen- 
tum, defined on continuous sets of configurations. 



CONTINUOUS HAMILTONIAN 

Rather than use Nd orthonormal quantum states to de- 
scribe a classical orbit with Nd informationally distinct 



configurations, it is sometimes convenient to use more. 
In the continuous-basis limit this yields a continuous- 
hamiltonian description. 



Oversampled dynamics 

Suppose that, starting with a classical finite-state dy- 
namics, we add M — 1 redundant intermediate-time states 
in the unit interval between each pair of consecutive 
integer-time states. Each orbit d of the corresponding 
quantum dynamics (generated by the hamiltonian Y\m) 
now visits MNd basis states rather than just the N d of 
the original dynamics (generated by Hi), and the state of 
the new Hm dynamics at a continuous moment of time 
becomes, from 

MN d -l 

\t,d,M}= ]T S(MN d ,k-Mt) \£,d,M) , (18) 

k=0 

where the basis state | jj , d, M) has been labeled by the 
time k/M when it is reached in an evolution starting 
from \0,d,M). Since this extended dynamics traverses 
distinct states at a rate um that is M times the original 
rate v, it has M times the average energy. As the num- 
ber of intermediate states added in a fixed time period 
goes to infinity, the hamiltonian Hm approaches a con- 
tinuous hamiltonian H^ and the average energy of the 
state \t,d,M) goes to infinity. 



Bandlimited basis 

By putting a bandlimit on the energy spectrum of the 
configuration basis states we can make the H m dynamics 
isomorphic to the original Hi dynamics, with the same 
average energy: a bandlimit on energy can correct for an 
oversampling of the underlying classical dynamics. 

The Fourier transform relationship ([2| between energy 
eigenstates and configurational basis states is left un- 
changed but we construct, in addition, a new set of band- 
limited configurations \n,d,M) N which are the Fourier 
transforms of the lowest Nd energy eigenstates of Hm, 
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(19) 

with n an integer. These states constitute an orthonor- 
mal basis for bandlimited superpositions of configura- 
tions. They have the same average energy as the con- 
figuration basis states of the Hi dynamics: the amount 
of time Td taken for one period of the orbit is being kept 
constant, and so from Q the first Nd energy eigenvalues 
mh/T d of Hm are the same as for Hi. 

The continuous time states \t, d, M) N that evolve from 
|0,d, M) N are given by (|10[) with n replaced by t. As in 
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^ they obey 
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\t,d,M) Nd = J2 S(N d ,n-t)\n,d,M) Nd 



(20) 



so the evolution of bandlimited states is isomorphic with 



that of \t, d). Moreover, from (19 1 with n t and ex- 
pressing \E: m,d,M) in terms of the MNj configura- 
tional basis states using ([2|, 
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\t,d,M) Nd = 



>M 



£ S(N d ,±-t)\±,d,M) 
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M-l N, 
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(21) 
(22) 



The bandlimited state is, at all times, an equally 
weighted superposition of M equivalent states, each of 
which corresponds to the extended classical configuration 
at time t represented in a different unit-time separated 
basis. Thus the correspondence of \t, d, M) N to classical 
configurations is the same as for \t, d). 

The state (21) is a sum over configurations separated 



in time by du = 1/M. If we normalize each configuration 
state to length \J~M instead of to length 1, this becomes 
delta-function normalization in the limit M — > oo and 



\t,d, 



du S(Nd, u — t) \u, d, oo ) 



From this and (151 



(t',d,oo\t,d,(X)) Nd = (i/,d\t,d) 



(23) 



(24) 



and so we can use the isomorphic \t, d) states to deter- 
mine amplitudes in the continuous configuration basis. 



PARTICLE MOTION 

A classical finite-state lattice dynamics is naturally de- 
scribed as a repeating sequence of invertible gate oper- 
ations 18J. In mapping this onto a quantum dynamics, 
the classical model can be implemented isomorphically 
as a sequence of local unitary operations. 

Fundamental physics is, however, normally described 
as particle dynamics. To make contact with this view- 
point we can recast finite-state lattice dynamics as parti- 
cle mechanics, following the motions of individual 1 's as 
if they were distinguishable particles. 



from lattice site to adjacent lattice site at a constant 
rate, with average speed v = 1. The motion is periodic 
in space, traversing N lattice sites in a distance L before 
repeating. At t = the particle is at x = 0. 

For this classical evolution, we can take the state of 
the system to be the integer position n of the 1 at inte- 
ger time n. In an isomorphic Hi quantum evolution, the 
distinct classical configurations become integer-position 
basis states \n). From (|9| we get a description of inter- 
mediate configurations in terms of integer-time ones, 
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where |a;) is the configuration obtained by evolving for 
a time t = x/v from the configuration |0). We identify 
the non-integer \x) with the non-integer positions of the 
continuously extended dynamics. 

In the quantum description of a classical particle at a 
non-integer position vt modulo L, there is some ampli- 
tude for the particle at more than one integer position. 



From (15) and using (24) we can interpret 



ip(x, t) = (x\vt) = S(N, x - vt) 



(26) 



to be the amplitude to find the particle at any continuous 
position x at time t, and compute the average momentum 
directly from i/;(x,f). 

Alternatively, we can instead start with an infinite- 
dimensional quantum hamiltonian that generates a con- 
tinuous shift in space in the +x direction at speed v: 



Hoc = vp } 



(27) 



with p = —ih-0- . 

" ox 

The direction of the shift is apparent from noting that 
HooV' = ihdip/dt implies ip(x,t) — tp(x — vt, 0) |29j . Now 
we can make this dynamics isomorphic to the H i discrete 
shift by bandlimiting the initial state so that the evolu- 
tion only traverses iV distinct states in the width L. Then 



from ( 23 1 the state corresponding to a classical particle at 
position vt in the position basis is ip{x, t) = S(N, x — vt), 



E 

P = — = 



(28) 



with energy E = hN/2T just as in Hi. From p7| , 

h 

2A ' 

where A = L/N. The state S(N, x—vt) achieves a general 
bound A > h/2p on the average separation of distinct 
states of a moving particle 

Since this description applies to any particle shifting 
uniformly in a lattice dynamics, (28) gives the corre- 
sponding momentum. Of course only lattice update op- 
erations that actually move a particle contribute to the 
shift-energy E = vp portion of its total energy [30] . 



Single particle 



Classical mechanics 



Consider a classical lattice dynamics in which a single 
particle, represented by a 1, hops in the +a; direction 



We can often consider a classical lattice-gas dynamics 
to be a discrete-time sampling of an idealized classical- 
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mechanical particle dynamics [5J [2T] that obeys Hamil- 
ton's equations, 



dH 

dq 3 



d Pj 
dt 



dH 

9pj 



dt 



(29) 



To make the lattice dynamics run faster by a factor k we 
reduce the interval between the discrete events, r — > t/k. 
From (29), this can be accomplished by letting H — > kH, 
which is exactly the energy scaling required by (J7|. 

We can't just rescale r arbitrarily while keeping the 
Pj's and qj's unchanged, however, because particle veloc- 
ities are limited by the speed of light. We can, instead, 
run the dynamics faster by putting the discrete events 
closer together in both time and space, leaving veloci- 
ties unchanged. If the distance between events A — > A/k, 
then the scale of the pj 's must be multiplied by k to get 
an overall scaling of H by k in ( 29 1 . This is exactly the 
momentum scaling required by 
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Indistinguishable particles 

Treating l's in a classical finite-state lattice dynamics 
as distinguishable particles — and keeping track of the dis- 
crete position and velocity of each 1 — dramatically over- 
represents the number of distinct states: all states with 
the same spatial pattern of l's and velocities correspond 
to a single state of the original lattice model. We can 
fix this over-representation in a quantum description of 
the distinghishable particle dynamics by merging equiva- 
lent states, adding them together to form new occupation 
number basis states, and using only these to describe the 
evolution. If we antisymmetrize each sum under parti- 
cle interchange, the new basis states will each have at 
most one 1 with a given position and velocity — we can 
symmetrize instead to allow more [19] . 

To describe a dynamics in which the number of ones 
changes with time, we can use creation and annihilation 
operators to add and remove particles from the state, 
while maintaining symmctrization. These field operators 
inherit fermionic or bosonic commutation rules from the 
symmetrization [51] . As we see from (251 (or from (20) 
for Hoc), a finite set of bandlimited basis states allows a 
particle to be added centered at any continuous position 
in space. In one dimension with one velocity, for example, 
the creation operator ty^(x) for any x is a superposition 
of the creation operators U^(n) for integer positions n, 



N-l 



n — x 



(30) 



Of course nothing essential is gained by using a con- 
tinuous space and time description, since a bandlimited 
continuous state is completely determined by its values 
at discrete positions and times. Similarly, nothing es- 
sential is gained by introducing fermionic field operators: 



there would be no need to maintain the antisymmetry of 
equivalent states if the original dynamics were described 
isomorphically in terms of local unitary operations [32] . 



UNCERTAINTY 



The particle described by (26) moves at a constant 



speed and is localized to a single position basis state of a 
finite-dimensional basis at all times (c/. |22|). This in no 
way conflicts with the uncertainty relations of quantum 
mechanics, which can be regarded as bounds on repre- 
senting information using limited bandwidth. 



Bandwidth bounds 

Constraints on time or position determine the mini- 
mum width of the energy or momentum eigenfrequency 
distribution needed to describe a state that meets the 
constraints. In the usual uncertainty bounds we also as- 
sociate a width with the time or position amplitude dis- 
tribution (23) . but in general other constraints on time 
or position can be used to determine a minimum width 
of energy or momentum eigenfrequencies. 

For example, suppose we have an exactly periodic evo- 
lution with period T. The state at time t can be written 
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Exact periodicity requires that each v n = E n /h be an in- 
teger multiple of 1/T. If this evolution passes through N 
mutually orthogonal states, then the superposition must 
involve at least N different |-E„)'s (since you can't con- 
struct N distinct states out of fewer than N distinct 
states). Moreover, there must also be at least N dis- 
tinct frequencies (since groups of |-E,i)'s with the same 
frequency act like a single eigenstate in the construc- 
tion). To have N distinct frequencies that are integer 
multiples of 1 /T, the bandwidth B (highest frequency in 
the superposition minus lowest) must obey 



B > 



N - 1 



(32) 



This is a version of the bandwidth-time theorem of com- 
munication theory [24 . If we let r = T/N be the aver- 
age time between distinct states, we see that this is also 
a version of the time-energy uncertainty relation, using 
B directly rather than some other measure of the width 
of the energy eigenfrequency distribution. The definition 
([3]) achieves this bound. 

Second-moment bounds 

In constructing uncertainty bounds, the standard de- 
viation of the eigenfrequency distribution is traditionally 
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chosen to measure its width. This choice reflects both 
familiarity from statistics and (for position and momen- 
tum) a simple connection between the commutation rela- 
tion and the standard-deviation bound [35J. This choice 
is often divergent, however, and so fails to provide a use- 
ful bound [53]. This is true in our case. 

Consider the bandlimited state (x\x) = S(N,x — x) 
centered at x. Limiting ourselves to spatial frequencies 
m/L with to ranging from to N — 1, this state has the 
least possible information about what the momentum is, 
since all momentum eigenstates in the allowed range have 
equal amplitude. Correspondingly we might expect the 
position to be as well-defined as possible, given the lim- 
ited bandwidth. It is clear from Figure[T]that the position 
localization of the probability distribution |i5(iV, x — x)\ 2 
is similar to that of a gaussian (dotted line). This is not 
apparent in the mean square position deviation, however, 



which can be estimated for large N using ( 11 ) as 



((x — a;) 2 ) ss / (x — x) 



N 



2 sin 7r(x — x) 

x — x I — 7T-, r^— ax 



n 2 (x — x) 2 



N 
2^2 



, (33) 



which diverges as N — > oo (i.e., on an infinitely wide 
space) [33] ■ Thus S(N,x — x), which is perfectly dis- 
tinct from a unit shift of itself, is not localized at all on 
the infinite line if we use the traditional second-moment 
measure of the width of the distribution. The unit-height 
gaussian, which looks so similar in the figure, has a mean 
square deviation of 1/27T. Other measures of the width 
have been proposed that avoid this disparity [23 . 



First-moment bounds 

For our purposes, a much better measure of the width 
of the eigenfrequency distribution is twice the average 
half-width: 2(9 — Here v is the average frequency 
[e.g., E/h) and Vq the lowest frequency used (e.g., Eg/h). 
In general [S], 



Uncertain states 

We have seen examples where a quantum hamiltonian 
describes a classical finite-state dynamics, but also makes 
extra distinctions not present in the original dynamics: 
A many particle hamiltonian that keeps track of which 
identical 1-bit is where. A continuous-shift hamiltonian 
that adds distinct states between the discrete time steps. 

We can eliminate over-representation and make the 
dynamics isomorphic to the original by adding together 
equivalent configurations with equal weight to construct 
truly distinct basis states. Starting from these, equiv- 
alent configurations will always have equal probability: 
equivalence is represented as uncertainty |34j . 

In the construction of the occupation number ba- 
sis states for identical 1-bits, a symmetrized or anti- 
symmetrized state represents equivalent states as being 
equally probable. In the case of over-representation of in- 
termediate states, constructing a basis without the high 
frequency information needed to represent intermediate 



details also merges equivalent states (21 1, making them 
equally probable. 

The continuous-hamiltonian representation of a dis- 
crete shift is an interesting limiting case of representing 
equivalence as uncertainty. A bandlimit with N distinct 
states yields (22) for finite M. For a state centered at 
t = x/v and M — > oo this becomes 

f 1 ( N ^ \ 

\x, oo) N = I du S(N, n + U — x) \n + u, oo) ] , 

J ° Vn=0 / 

(35) 



which is a uniform superposition of all the equivalent 
ways to represent a classical particle at position x if only 
N equally-spaced positions are distinct. 

The tradeoff between bandwidth and minimum sep- 
aration in space determines the minimum uncertainty 
volume of phase space needed to represent each distinct 
state [35], and this is achieved by \x,oo) N . 



2(9 - vq) > B a 



(34) 



DISCUSSION 



where B m - m is the minimum bandwidth compatible with 
the temporal or spatial constraints on the system. 

For example, if r mm is the minimum separation in time 
between two mutually orthogonal states in the evolution, 
then the minimum bandwidth needed is B min — l/2r m i n : 
there must be at least two distinct frequencies and they 
must be separated by at least half of 1/T mm . 

The B m in = l/2r m i n bound (34) is only achieved by 
the energy ^ for N = 2. For TV > 2, the energy @ 
is about twice as great as allowed by this bound. There 



is, however, the additional bandwidth constraint (32) re- 



quired to have N distinct states in period T. The energy 
^ achieves (34) with this constraint. 



Classical finite-state dynamics that are invertible can 
be mapped isomorphically onto the discrete time behav- 
ior of finite-energy quantum dynamics. A quantum evo- 
lution mapping an infinite number of distinct states into 
a finite time period would have an infinite average energy. 

Quantum-classical isomorphism challenges conven- 
tional wisdom about essential differences between quan- 
tum and classical systems: identical particles, ampli- 
tudes, frequencies, complementarity and uncertainty all 
play essential roles in describing and analyzing classical 
finite-state dynamics using continuous language. 

Quantum-classical models also shed light on the foun- 
dations of classical mechanics. They provide a quantum 



substrate where interesting classical behavior arises with- 
out approximation or decoherence. Physically meaning- 
ful energy and momentum scales are defined directly by 
the separation of classical events in time and space. 

Finally, quantum-classical isomorphism may be helpful 
in understanding and teaching quantum mechanics. Just 
as it is useful to study classical information and classical 
computation as a preliminary to studying their quantum 
counterparts, it seems useful to study other aspects of 
the machinery and concepts of quantum mechanics in a 
simplified classical setting. 
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This is equivalent to phase and number [8] Eqn. 41] , with 
9\n,d) = (2%n/N d ) \n,d) and n\E: m,d) = m\E: m,d). 
In this example the local additivity of the energy derives 
from the locality of the sequence of update operations; 
H doesn't need to be expressible in a manifestly local 
fashion (cf. [10] [11]). 

This idea is the basis of practical lattice gas models of 
classical hydrodynamics [12] and of quantum fluids [17] 
(both with continuous isotropy in the macroscopic limit). 
For motion in the —x direction we would use Hoo = — vp 
instead. To represent the direction explicitly in the state 
we would multiply it by [+) or |— ) and let Hoo = o z vp. 
To make updates local, we can use a partitioning dynam- 
ics, in which each particle motion involves a single update 
that changes both old and new positions at once [4l I18|. 
For example, if particle labels are generated sequentially 
as particles are created, then interchanging the order 
in which two particles are added to an antisymmetrized 
state is equivalent to interchanging their particle labels, 
and so creation operators must anticommute [20] . 
This should apply equally to quantum lattice gas simu- 
lations of non-classical systems. 

If x is near the middle of the periodic space, then the 
wavefunction goes to at the boundary as N — > oo and 
so the usual uncertainty relations apply 8, Eqn. 22]. 
Ignorance of differences between equivalent states doesn't 
count toward entropy, which is one reason quantum prob- 
abilities must be kept separate from ordinary ones [4]. 
In a periodic space of length L, momentum eigenfrequen- 
cies must be integer multiples of 1/L. Thus to represent 
N distinct states a bandwidth B > (N — 1)/L is needed, 
and so the frequency-space volume per distinct state is 
BL/N > (N-l)/N (uncertainty tradeoff for B vs. L/N). 



